Genomic analysis supports Cape Lion population connectivity prior to colonial eradication and extinction.

Cape lions (Panthera leo melanochaitus) formerly ranged throughout the grassland plains of the "Cape Flats" in what is today known as the Western Cape Province, South Africa. Cape lions were likely eradicated because of overhunting and habitat loss after European colonization. European naturalists originally described Cape lions as "black-maned lions" and claimed that they were phenotypically distinct. However, other depictions and historical descriptions of lions from the Cape report mixed or light coloration and without black or extensively developed manes. These findings suggest that, rather than forming a distinct population, Cape lions may have had phenotypic and genotypic variation similar to other African lions. Here we investigate Cape lion genome characteristics, population dynamics, and genetic distinctiveness prior to their extinction. We generated genomic data from 2 historic Cape lions to compare to 118 existing high-coverage mitogenomes, and low-coverage nuclear genomes of 53 lions from 13 African countries. We show that, before their eradication, lions from the Cape Flats had diverse mitogenomes and nuclear genomes that clustered with lions from both southern and eastern Africa. Cape lions had high genome-wide heterozygosity and low inbreeding coefficients, indicating that populations in the Cape Flats went extinct so rapidly that genomic effects associated with long-term small population size and isolation were not detectable. Our findings do not support the characterization of Cape lions as phylogeographically distinct, as originally put forth by some European naturalists, and illustrates how alternative knowledge systems, for example, Indigenous perspectives, could potentially further inform interpretations of species histories.

Lion Localizer: a software tool for inferring the provenance of lions (Panthera leo) using mitochondrial DNA.

The illegal poaching of lions for their body parts poses a severe threat to lion populations across Africa. Poaching accounts for 35% of all human-caused lion deaths, with 51% attributed to retaliatory killings following livestock predation. In nearly half of the retaliatory killings, lion body parts are removed, suggesting that high demand for lion body parts may fuel killings attributed to human-lion conflict. Trafficked items are often confiscated in transit or destination countries far from their country of origin. DNA from lion parts may in some cases be the only available means for examining their geographic origins. In this paper, we present the Lion Localizer, a full-stack software tool that houses a comprehensive database of lion mitochondrial DNA (mtDNA) sequences sourced from previously published studies. The database covers 146 localities from across the African continent and India, providing information on the potential provenance of seized lion body parts. Lion mtDNA sequences of 350 bp or 1140 bp corresponding to the cytochrome b region can be generated from lion products and queried against the Lion Localizer database. Using the query sequence, the Lion Localizer generates a listing of exact or partial matches, which are displayed on an interactive map of Africa. This allows for the rapid identification of potential regions and localities where lions have been or are presently being targeted by poachers. By examining the potential provenance of lion samples, the Lion Localizer serves as a valuable resource in the fight against lion poaching. The software is available at https : //lionlocalizer.org.

Swine global genomic resources: insights into wild and domesticated populations.

Suids, both domesticated and wild, are found on all continents except for Antarctica and provide valuable food resources for humans in addition to serving as important models for biomedical research. Continuing advances in genome sequencing have allowed researchers to compare the genomes from diverse populations of suids helping to clarify their evolution and dispersal. Further analysis of these samples may provide clues to improve disease resistance/resilience and productivity in domestic suids as well as better ways of classifying and conserving genetic diversity within wild and captive suids. Collecting samples from diverse populations of suids is resource intensive and may negatively impact endangered populations. Here we catalog extensive tissue and DNA samples from suids in collections in both Europe and North America. We include samples that have previously been used for whole genome sequencing, targeted DNA sequencing, RNA sequencing, and reduced representation bisulfite sequencing (RRBS). This work provides an important centralized resource for researchers who wish to access published databases.

Evolution: Untangling the woolly mammoth.

Twenty-two woolly mammoth genomes have been compared to those of living elephants, identifying genes under strong evolutionary pressure in mammoths, including genes associated with curly, wiry, thick, bushy, coarse, uncombable and (of course) woolly hair.

The Impact of Variation in the Toll-like Receptor 3 Gene on Epizootic Hemorrhagic Disease in Illinois Wild White-Tailed Deer (Odocoileus virginianus).

Epizootic hemorrhagic disease (EHD) leads to high mortality in white-tailed deer (Odocoileus virginianus) and is caused by a double-stranded RNA (dsRNA) virus. Toll-like receptor 3 (TLR3) plays a role in host immune detection and response to dsRNA viruses. We, therefore, examined the role of genetic variation within the TLR3 gene in EHD among 84 Illinois wild white-tailed deer (26 EHD-positive deer and 58 EHD-negative controls). The entire coding region of the TLR3 gene was sequenced: 2715 base pairs encoding 904 amino acids. We identified 85 haplotypes with 77 single nucleotide polymorphisms (SNPs), of which 45 were synonymous mutations and 32 were non-synonymous. Two non-synonymous SNPs differed significantly in frequency between EHD-positive and EHD-negative deer. In the EHD-positive deer, phenylalanine was relatively less likely to be encoded at codon positions 59 and 116, whereas leucine and serine (respectively) were detected less frequently in EHD-negative deer. Both amino acid substitutions were predicted to impact protein structure or function. Understanding associations between TLR3 polymorphisms and EHD provides insights into the role of host genetics in outbreaks of EHD in deer, which may allow wildlife agencies to better understand the severity of outbreaks.

Numt Parser: Automated identification and removal of nuclear mitochondrial pseudogenes (numts) for accurate mitochondrial genome reconstruction in Panthera.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from 2 ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to 2 other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

County-wide assessments of Illinois white-tailed deer (Odocoileus virginianus) prion protein gene variation using improved primers and potential implications for management.

Chronic wasting disease (CWD) is a fatal, highly infectious prion disease that affects captive and wild cervids. Chronic wasting disease is the only known transmissible spongiform encephalopathy affecting free-ranging wildlife. In CWD-positive deer, some haplotypes of the prion protein gene PRNP are detected at lower frequencies as compared to CWD-negative deer, as are some variants of the prion protein PrP. Here, we examined wild, hunter-harvested CWD-negative white-tailed deer (Odocoileus virginianus) to determine whether there were geographical or temporal differences in the PRNP haplotypes, PRNP diplotypes, PrP proteoforms, and in the proportion of deer with at least one protective haplotype. We sampled 96-100 hunter-harvested deer per county at two time points in the Illinois counties of Jo Daviess, LaSalle, and Winnebago, chosen based on their geographic locations and known occurrence of CWD. The entire coding region of PRNP was sequenced, with haplotypes, diplotypes, and PrP proteoforms inferred. Across time, in Winnebago there was a significant increase in PrP proteoform F (p = 0.034), which is associated with a lower vulnerability to CWD. In every county, there was an increase over time in the frequency of deer carrying at least one protective haplotype to CWD, with a significant increase (p = 0.02) in the Jo Daviess County CWD infected region. We also found that primer combination was important as there was an 18.7% difference in the number of the deer identified as homozygous depending on primer usage. Current Illinois state management practices continue to remove CWD infected deer from locally infected areas helping to keep CWD prevalence low. Nonetheless, continued research on spatial and temporal changes in PRNP haplotypes, PrP proteoforms, and levels of deer vulnerability among Illinois deer will be important for the management of CWD within the state of Illinois and beyond.

A De Novo Chromosome-Level Genome Assembly of the White-Tailed Deer, Odocoileus Virginianus.

Cervids are distinguished by the shedding and regrowth of antlers. Furthermore, they provide insights into prion and other diseases. Genomic resources can facilitate studies of the genetic underpinnings of deer phenotypes, behavior, and disease resistance. Widely distributed in North America, the white-tailed deer (Odocoileus virginianus) has recreational, commercial, and food source value for many households. We present a genome generated using DNA from a single Illinois white-tailed sequenced on the PacBio Sequel II platform and assembled using Wtdbg2. Omni-C chromatin conformation capture sequencing was used to scaffold the genome contigs. The final assembly was 2.42 Gb, consisting of 508 scaffolds with a contig N50 of 21.7 Mb, a scaffold N50 of 52.4 Mb, and a BUSCO complete score of 93.1%. Thirty-six chromosome pseudomolecules comprised 93% of the entire sequenced genome length. A total of 20 651 predicted genes using the BRAKER pipeline were validated using InterProScan. Chromosome length assembly sequences were aligned to the genomes of related species to reveal corresponding chromosomes.

Spatial epidemiology of hemorrhagic disease in Illinois wild white-tailed deer.

Epizootic hemorrhagic disease (EHD) and bluetongue (BT) are vector-borne viral diseases that affect wild and domestic ruminants. Clinical signs of EHD and BT are similar; thus, the syndrome is referred to as hemorrhagic disease (HD). Syndromic surveillance and virus detection in North America reveal a northern expansion of HD. High mortalities at northern latitudes suggest recent incursions of HD viruses into northern geographic areas. We evaluated the occurrence of HD in wild Illinois white-tailed deer from 1982 to 2019. Our retrospective space-time analysis identified high-rate clusters of HD cases from 2006 to 2019. The pattern of northward expansion indicates changes in virus-host-vector interactions. Serological evidence from harvested deer revealed prior infection with BTV. However, BTV was not detected from virus isolation in dead deer sampled during outbreaks. Our findings suggest the value of capturing the precise geographic location of outbreaks, the importance of virus isolation to confirm the cause of an outbreak, and the importance of expanding HD surveillance to hunter-harvested wild white-tailed deer. Similarly, it assists in predicting future outbreaks, allowing for targeted disease and vector surveillance, helping wildlife agencies communicate with the public the cause of mortality events and viral hemorrhagic disease outcomes at local and regional scales.

Combining methods for non-invasive fecal DNA enables whole genome and metagenomic analyses in wildlife biology.

Non-invasive biological samples benefit studies that investigate rare, elusive, endangered, or dangerous species. Integrating genomic techniques that use non-invasive biological sampling with advances in computational approaches can benefit and inform wildlife conservation and management. Here, we used non-invasive fecal DNA samples to generate low- to medium-coverage genomes (e.g., >90% of the complete nuclear genome at six X-fold coverage) and metagenomic sequences, combining widely available and accessible DNA collection cards with commonly used DNA extraction and library building approaches. DNA preservation cards are easy to transport and can be stored non-refrigerated, avoiding cumbersome or costly sample methods. The genomic library construction and shotgun sequencing approach did not require enrichment or targeted DNA amplification. The utility and potential of the data generated was demonstrated through genome scale and metagenomic analyses of zoo and free-ranging African savanna elephants (Loxodonta africana). Fecal samples collected from free-ranging individuals contained an average of 12.41% (5.54-21.65%) endogenous elephant DNA. Clustering of these elephants with others from the same geographic region was demonstrated by a principal component analysis of genetic variation using nuclear genome-wide SNPs. Metagenomic analyses identified taxa that included Loxodonta, green plants, fungi, arthropods, bacteria, viruses and archaea, showcasing the utility of this approach for addressing complementary questions based on host-associated DNA, e.g., pathogen and parasite identification. The molecular and bioinformatic analyses presented here contributes towards the expansion and application of genomic techniques to conservation science and practice.

Evaluating the ability of a locally focused culling program in removing chronic wasting disease infected free-ranging white-tailed deer in Illinois, USA, 2003-2020.

In northern Illinois, chronic wasting disease (CWD) was first identified in free-ranging white-tailed deer (Odocoileus virginianus; hereafter referred to as "deer") in 2002. To reduce CWD transmission rates in Illinois, wildlife biologists have conducted locally focussed culling of deer since 2003 in areas where CWD has been detected. We used retrospective spatial, temporal and space-time scan statistical models to identify areas and periods where culling removed higher than expected numbers of CWD-positive deer. We included 490 Public Land Survey "sections" (∼2.59 km2 ) from 15 northern Illinois counties in which at least one deer tested positive for CWD between 2003 and 2020. A negative binomial regression model compared the proportion of CWD positive cases removed from sections with at least one CWD case detected in the previous years, "local area 1 (L1)," to the proportion of CWD cases in adjacent sections-L2, L3, and L4-designated by their increasing distance from L1. Of the 14,661 deer removed and tested via culling, 325 (2.22 %) were CWD-positive. A single temporal CWD cluster occurred in 2020. Three spatial clusters were identified, with a primary cluster located at the border of Boone and Winnebago counties. Four space-time clusters were identified with a primary cluster in the northern portion of the study area from 2003 to 2005 that overlapped with the spatial cluster. The proportion of CWD cases removed from L1 (3.92, 95% CI, 2.56-6.01) and L2 (2.32, 95% CI, 1.50-3.59) were significantly higher compared to L3. Focussing culling efforts on accessible properties closest to L1 areas results in more CWD-infected deer being removed, which highlights the value of collaborations among landowners, hunters, and wildlife management agencies to control CWD. Continuous evaluation and updating of the culling and surveillance programs are essential to mitigate the health burden of CWD on deer populations in Illinois.

Bluetongue and Epizootic Hemorrhagic Disease in the United States of America at the Wildlife-Livestock Interface.

Bluetongue (BT) and epizootic hemorrhagic disease (EHD) cases have increased worldwide, causing significant economic loss to ruminant livestock production and detrimental effects to susceptible wildlife populations. In recent decades, hemorrhagic disease cases have been reported over expanding geographic areas in the United States. Effective BT and EHD prevention and control strategies for livestock and monitoring of these diseases in wildlife populations depend on an accurate understanding of the distribution of BT and EHD viruses in domestic and wild ruminants and their vectors, the Culicoides biting midges that transmit them. However, national maps showing the distribution of BT and EHD viruses and the presence of Culicoides vectors are incomplete or not available at all. Thus, efforts to accurately describe the potential risk of these viruses on ruminant populations are obstructed by the lack of systematic and routine surveillance of their hosts and vectors. In this review, we: (1) outline animal health impacts of BT and EHD in the USA; (2) describe current knowledge of the distribution and abundance of BT and EHD and their vectors in the USA; and (3) highlight the importance of disease (BT and EHD) and vector surveillance for ruminant populations.

Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis).

OBJECTIVE: The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling. RESULTS: Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; HO = 0.30). These were sufficient to distinguish among individuals (PID < 0.0001), and to distinguish among siblings (PID(sib) < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.

Retroviral integrations contribute to elevated host cancer rates during germline invasion.

Repeated retroviral infections of vertebrate germlines have made endogenous retroviruses ubiquitous features of mammalian genomes. However, millions of years of evolution obscure many of the immediate repercussions of retroviral endogenisation on host health. Here we examine retroviral endogenisation during its earliest stages in the koala (Phascolarctos cinereus), a species undergoing germline invasion by koala retrovirus (KoRV) and affected by high cancer prevalence. We characterise KoRV integration sites (IS) in tumour and healthy tissues from 10 koalas, detecting 1002 unique IS, with hotspots of integration occurring in the vicinity of known cancer genes. We find that tumours accumulate novel IS, with proximate genes over-represented for cancer associations. We detect dysregulation of genes containing IS and identify a highly-expressed transduced oncogene. Our data provide insights into the tremendous mutational load suffered by the host during active retroviral germline invasion, a process repeatedly experienced and overcome during the evolution of vertebrate lineages.

Spatial analysis of chronic wasting disease in free-ranging white-tailed deer (Odocoileus virginianus) in Illinois, 2008-2019.

Understanding the geographic distribution and clustering of chronic wasting disease (CWD) among free-ranging white-tailed deer (Odocoileus virginianus) populations can inform disease management practices. We used a retrospective analysis of surveillance data to evaluate CWD's spatial and temporal dynamics within 16 CWD-infected northern Illinois counties. Of 42,541 deer samples collected and tested for CWD from recreational hunter harvest between 2008 and 2019, we recorded 359 (0.84%) CWD-positive samples. We observed variability in CWD cases over time and space. By county, the median CWD-positive proportion was 0.84%, varying from a minimum of 0.14% in McHenry County to a maximum of 6.28% in Boone County. Across years, there were differences among CWD-positive proportions with a median of 0.90%, ranging from a minimum of 0.27% in 2012 to a maximum of 1.60% in 2019. We used a retrospective discrete Poisson scan statistic model to evaluate the space-time clustering of CWD-positive deer. We identified a statistically significant (p < .001) primary cluster C1 (area = 23.59 km2 ; RR = 10.48), occurring from 2010 to 2015 in the north-central part of the study area, and a secondary cluster C2, occurring from 2014 to 2019 (area = 9.27 km2 ; RR = 3.88) in the north-west of the study area. Detected CWD-positive space-time clusters suggest that the risk of CWD is not random. Space-time clusters of CWD can be used to evaluate the effectiveness of the Illinois CWD management programme. The area surrounding the older C1 cluster has undergone longer and more intense CWD management compared with C2. Currently, the older C1 cluster is no longer as high risk compared with the newer cluster C2, suggesting that management efforts in C2 should be increased. However, all CWD clusters should be targeted with surveillance, prevention and management programmes, including reducing deer densities to limit further spread of CWD.

Sourcing Elephant Ivory from a Sixteenth-Century Portuguese Shipwreck.

The oldest known shipwreck in southern Africa was found in Namibia in 2008.1-4 Forty tons of cargo, including gold and silver coins, helped identify the ship as the Bom Jesus, a Portuguese nau (trading vessel) lost in 1533 while headed to India.4-6 The cargo included >100 elephant tusks,7 which we examined using paleogenomic and stable isotope analyses. Nuclear DNA identified the ivory source as African forest (Loxodonta cyclotis) rather than savanna (Loxodonta africana) elephants. Mitochondrial sequences traced them to West and not Central Africa and from ≥17 herds with distinct haplotypes. Four of the haplotypes are known from modern populations; others were potentially lost to subsequent hunting of elephants for ivory. Stable isotope analyses (δ13C and δ15N) indicated that the elephants were not from deep rainforests but from savanna and mixed habitats. Such habitats surround the Guinean forest block of West Africa8 and accord with the locations of major historic Portuguese trading ports.9,10 West African forest elephants currently range into savanna habitats;11-13 our findings suggest that this was not consequent to regional decimation of savanna elephants for their ivory in the 19th and 20th centuries. During the time of the Bom Jesus, ivory was a central driver in the formation of maritime trading systems connecting Europe, Africa, and Asia. Our integration of paleogenomic, archeological, and historical methods to analyze the Bom Jesus ivory provides a framework for examining vast collections of archaeological ivories around the world, in shipwrecks and other contexts.

Evolution: Untangling the Woolly Rhino's Extinction.

The woolly rhinoceros was a charismatic inhabitant of the frigid steppes of Pleistocene Eurasia. Now, the genome of an 18,500-year-old woolly rhino has been sequenced. It points to a thriving population less than 5000 years before the species disappeared.

Prion Protein Gene (PRNP) Sequences Suggest Differing Vulnerability to Chronic Wasting Disease for Florida Key Deer (Odocoileus virginianus clavium) and Columbian White-Tailed Deer (O. v. leucurus).

Chronic wasting disease (CWD) is a fatal, highly transmissible spongiform encephalopathy caused by an infectious prion protein. CWD is spreading across North American cervids. Studies of the prion protein gene (PRNP) in white-tailed deer (WTD; Odocoileus virginianus) have identified non-synonymous substitutions associated with reduced CWD frequency. Because CWD is spreading rapidly geographically, it may impact cervids of conservation concern. Here, we examined the genetic vulnerability to CWD of 2 subspecies of WTD: the endangered Florida Key deer (O. v. clavium) and the threatened Columbian WTD (O. v. leucurus). In Key deer (n = 48), we identified 3 haplotypes formed by 5 polymorphisms, of which 2 were non-synonymous. The polymorphism c.574G>A, unique to Key deer (29 of 96 chromosomes), encodes a non-synonymous substitution from valine to isoleucine at codon 192. In 91 of 96 chromosomes, Key deer carried c.286G>A (G96S), previously associated with substantially reduced susceptibility to CWD. Key deer may be less genetically susceptible to CWD than many mainland WTD populations. In Columbian WTD (n = 13), 2 haplotypes separated by one synonymous substitution (c.438C>T) were identified. All of the Columbian WTD carried alleles that in other mainland populations are associated with relatively high susceptibility to CWD. While larger sampling is needed, future management plans should consider that Columbian WTD are likely to be genetically more vulnerable to CWD than many other WTD populations. Finally, we suggest that genetic vulnerability to CWD be assessed by sequencing PRNP across other endangered cervids, both wild and in captive breeding facilities.

Association of chronic wasting disease susceptibility with prion protein variation in white-tailed deer (Odocoileus virginianus).

Chronic wasting disease (CWD) is caused by prions, infectious proteinaceous particles, PrPCWD. We sequenced the PRNP gene of 2,899 white-tailed deer (WTD) from Illinois and southern Wisconsin, finding 38 haplotypes. Haplotypes A, B, D, E, G and 9 others encoded Q95G96S100N103A123Q226, designated 'PrP variant A.' Haplotype C and 4 other haplotypes encoded PrP 'variant C' (Q95S96S100N103A123Q226). Haplotype F and two other haplotypes encoded PrP 'variant F' (H95G96S100N103A123Q226). The association of CWD with encoded PrP variants was examined in 2,537 tested WTD from counties with CWD. Relative to PrP variant A, CWD susceptibility was lower in deer with PrP variant C (OR = 0.26, p < 0.001), and even lower in deer with PrP variant F (OR = 0.10, p < 0.0001). Susceptibility to CWD was highest in deer with both chromosomes encoding PrP variant A, lower with one copy encoding PrP variant A (OR = 0.25, p < 0.0001) and lowest in deer without PrP variant A (OR = 0.07, p < 0.0001). There appeared to be incomplete dominance for haplotypes encoding PrP variant C in reducing CWD susceptibility. Deer with both chromosomes encoding PrP variant F (FF) or one encoding PrP variant C and the other F (CF) were all CWD negative. Our results suggest that an increased population frequency of PrP variants C or F and a reduced frequency of PrP variant A may reduce the risk of CWD infection. Understanding the population and geographic distribution of PRNP polymorphisms may be a useful tool in CWD management.